A closer look at arrested spinodal decomposition in protein solutions 
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Concentrated aqueous solutions of the protein lysozyme undergo a liquid-solid transition upon a 
temperature quench into the unstable spinodal region below a characteristic arrest temperature of 
Tf = 15 °C. We use video microscopy and ultra small angle light scattering in order to investigate 
the arrested structures as a function of initial concentration, quench temperature and rate of the 
temperature quench. We find that the solid-like samples show all the features of a bicontinuous 
network that is formed through a spinodal decomposition process. We determine the correlation 
length £ and demonstrate that £ exhibits a temperature dependence that closely follows critical 
scaling expected for density fluctuations during the early stages of spinodal decomposition. These 
findings are in agreement with an arrest scenario based on a state diagram where the arrest or gel 
line extends far into the unstable region below the spinodal line. Arrest then occurs when during 
the early stage of spinodal decomposition the volume fraction (j>2 of the dense phase intersects the 
dynamical arrest threshold <f>2,Giass, upon which phase separation gets pinned into a space spanning 
gel network with a characteristic length £. 

PACS numbers: pacs 



I. INTRODUCTION 

Gelation of colloidal particles has continued to attract 
considerable attention during the last few years. While it 
is ubiquitous and used in numerous applications in ma- 
terial science, food science and the industry of personal 
care products 1 , it still remains a quite controversial 
topic and only partially understood. Under the general 
head line of "Dynamical arrest" or "Jamming" , attempts 
have been made to arrive at a much more fundamen- 
tal understanding of various fluid-solid transitions in col- 
loidal suspensions that encompass seemingly unrelated 
phenomena such as the glass transition observed for hard 
spheres, and gelation in suspensions of strongly attrac- 
tive and thus irreversibly aggregating particles. Particu- 
lar attention has been given to particles interacting via a 
short-ranged attractive potential, and the role of the key 
parameters volume fraction </>, interparticle interaction 
strength U a and range of the potential A in determining 
the state diagram and in particular the gelation line has 
been addressed [HOI!]. 

Colloids interacting via short-ranged attractions ex- 
hibit a number of interesting phenomena such as a re- 
entrant glass transition at high values of <fi and and a 
competition between phase separation and dynamical ar- 
rest at intermediate values of U a /kT and 4> that can lead 
to the phenomenon of an arrested spinodal decomposi- 
tion resulting in the formation of a solild-like network 
[BEEmHEIlinidl- Moreover, a recent study of Lu 
et al. provided evidence that gelation for short range at- 
tractive particles would occur as a consequence of an ini- 
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tial equilibrium liquid-liquid phase separation that would 
then drive the formation of a space-spanning cluster that 
dynamically arrest to form a gel [7j. However, while the 
phenomenon of an arrested spinodal decomposition as 
a possible mechanism for gelation in attractive colloidal 
suspensions appears to be generally accepted, there exists 
a controversy as to the location of the arrest line [Sj [7]. 

The presence of an arrested spinodal decomposition 
has for example been demonstrated for concentrated so- 
lutions of the model protein lysozyme, which is known to 
exhibit a phase diagram that closely matches the predic- 
tions for colloids with short-ranged attractions [5]. Here 
the mechanism underlying the liquid-solid transition en- 
countered when quenching a protein sample deeply below 
the coexistence curve for liquid-liquid phase separation 
has been interpreted based on a state diagram shown in 
Figjl] Here the arrest or gelation line extends deep into 
the unstable region, and arrest occurs as the spinodal 
decomposition process leads to a bi-continuous structure 
which gets " pinned" into a rigid self-supporting network 
when the concentration of the high density regions crosses 
the arrest line and subsequently undergo dynamical ar- 
rest. The findings described by Cardinaux et al. [5] are 
however in sharp contrast to one of the central conclu- 
sions drawn in ref. [7], where a universal phase dia- 
gram for colloids with short-range attraction had been 
constructed that suggested that the gelation line coin- 
cides with the phase separation boundary in the Baxter 
model. These authors then concluded that the origin of 
dynamic arrest came from the dense phase undergoing 
an attractive glass transition at cf> g ~ 0.55, and claimed 
that the attractive glass line would thus not extend into 
the phase separation region but instead follow its high 
density boundary. 

It is clear that this apparent discrepancy warrants fur- 
ther investigations. Moreover, given the fact that the lo- 
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FIG. 1: (a) State diagram of aqueous lysozyme solutions 
showing regions of crystal formation (open circles: experi- 
mentally determined phase boundary), and the metastable 
liquid-liquid phase separation below the experimentally de- 
terimined coexistence curve (open squares). The region below 
the coexistence curve can be separated into three areas differ- 
ing in their kinetic behavior: A region of complete demixing 
(I), gel formation through an arrested spinodal decomposi- 
tion (II), and a homogeneous attractive glass (III). Full sym- 
bols stand for the results of the centrifugation experiments 
described in [5] . (b) Temporal evolution of the sample tem- 
perature and the concomitant values of the storage and loss 
moduli G' and G", respectively at a frequency of xx Hz dur- 
ing a temperature quench as indicated by the arrow in Fig. 1 
a. (c) Schematic description of the model used to understand 
the behavior observed in Fig. 1 a and b. It schematically 
describes the evolution of density fluctuations as a function 
of time (ti < t2 < ts) during the early stage of spinodal de- 
composition. The initial and equilibrium volume fractions of 
the dense and dilute phases are denoted 4>o, <i>i and cj>\. <j>2, eq 
and <j>i, eq are the volume fraction of the right and left branch 
of the coexistence curve. A is the wavelength of the density 
fluctuations. At a certain moment tz the volume fraction (f>2 
of the dense phase intersects the dynamical arrest threshold 
4>2,Giass upon which phase separation gets pinned into a space 
spanning gel network with a characteristic length A 



cation of the glass line in the experiments with lysozyme 
had been determined indirectly by exploiting the unusual 
rheological properties of the arrested samples, these find- 
ings need further support. We thus now aim at investi- 
gating the sequence of events that finally lead to an ar- 
rested spinodal decomposition in concentrated solutions 
of lysozyme, and at obtaining an improved understand- 
ing of the link between quench depth, initial concentra- 
tion and the resulting structural and mechanical prop- 
erties. The main focus of the present paper is on the 
determination the characteristic length scale of the struc- 
tures formed in the arrested spinodal decomposition as a 
function of <p and final quench temperature Tf, and the 
attempt to explain the results based on known univer- 
sal scaling behavior of spinodal decomposition combined 
with an arrest scenario based on the phase diagram pos- 
tulated by Cardinaux et al. and reproduced again in Fig. 
[I] In the first part we present the structure of lysozyme 



samples quenched into the region between the spinodal 
line and the tie line at 15°C. We demonstrate that these 
samples exhibit indeed the classical features of spinodal 
decomposition. We are able to quantitatively follow the 
temporal evolution of the transient bicountinous network 
during the so-called intermediate regime, determine the 
characteristic length £ of the network from a combina- 
tion of video microscopy and ultra-small angle light scat- 
tering, and verify that it follows the expected universal 
scaling relationship where £ evolves with time t as t 1 / 3 . 

In the second part we then explore the structure of 
samples quenched in the spinodal region below the arrest 
tie line at 15°C. While we also observe the initial forma- 
tion of a bicountinous network, the structure becomes 
rapidly time-independent. We show that £ now depends 
mainly on the quench rate and the quench depth, but 
not on the volume fraction of lysozyme, and that the 
data are consistent with linear Cahn-Hillard theory, thus 
supporting a mechanism where the sample arrest in the 
early stage of spinodal decomposition consistent with the 
position of the dynamical arrest line as determined pre- 
viously [5]. 



II. MATERIAL AND METHODS 

We use hen egg white lysozyme (Fluka, L7651) dis- 
persed in an aqueous buffer (20 mM Hepes) containing 
0.5 M sodium chloride. Details of the sample prepara- 
tion procedures can be found elsewhere EI: . Initially 
a suspension at (/) ~ 0.22 is prepared in pure buffer with- 
out added salt, and its pH is adjusted to 7.8 ± 0.1 with 
sodium hydroxide j5J H3] . We then dilute it with a NaCl- 
containing buffer at pH= 7.8 to a final NaCl concen- 
tration of 0.5 M. Particular care is taken to avoid partial 
phase separation upon mixing by pre-heating both buffer 
and stock solution well above the coexistence curve for 
liquid- liquid phase separation (cf. Fig. 1). This proce- 
dure results in completely transparent samples at room 
temperature with volume fractions ranging from 4> = 0.01 
to 0.19, where <f> is obtained from the protein concentra- 
tions c measured by UV absorption spectroscopy using 
cf> = c/p, where p — 1.351 g/cm 3 is the protein partial 
specific density, [13] . 

We use a combination of ultra small angle light scatter- 
ing (USALS) and video microscopy to study the temporal 
evolution of the bicontinuous network structure during 
spinodal decomposition of the lysozyme solutions. We 
investigate two types of temperature quenches differing 
only by their cooling rate. For both types of quenches, a 
fresh lysozyme dispersion is filled in a cuvette at a tem- 
perature around 25° C, and the sample is then quenched 
below the spinodal line to a final temperature Tf. We 
use either rectangular scattering cells with a short path 
length between 10 /im and 20 pm to avoid multiple scat- 
tering for the USAL or capillary tubes (width 1 mm, 
depth 50 /im) for video microscopy. In the 'fast' quench 
the sample is pre-quenched for a minute in an ethanol 
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bath at Tf. The sample is then quickly transferred to 
the thermostated cuvette holder of the experiment, also 
at Tf . This procedure allows us to obtain the fastest ex- 
perimental quenches possible with the cell size required 
by the experimental set-up. We estimate the time it takes 
for the sample to go from 25° C to Tf to be less than 30 
s. 'Slow' quenches result from a simple cooling of the 
sample from 25° C to Tf as the sample is inserted in the 
pre-thermostated sample holder of the experiment. It 
then takes about 100 s to thermalize the sample to Tf. 

The USALS setup covers a range of scattering vectors 
0.1<<? <2 /jm" 1 . It is described in detail elsewhere [15]. 
Video microscopy is performed with a Leica DM- IRB 
in phase contrast mode. To minimize heating during 
observation, we illuminate the sample only during the 
time required to take pictures. Moreover we attenuate 
the heating effect of the UV and IR part of the illumi- 
nation spectrum by placing a beaker of water between 
the light source and the sample. The microscope focus 
is set to image the mid- plane of the sample to avoid 
the influence from wetting effects at the interfaces be- 
tween the lysozyme solution and the capillary, and the 
objective aperture is opened to its maximum diameter 
to minimize the depth of the focal plan to a few fim. We 
carefully monitor the sample for crystallization and dis- 
regard all samples where crystals form during the course 
of the experiment. 



III. 



'CLASSICAL' PHASE SEPARATION 



For shallow quenches below the spinodal curve but 
above the arrest tie line at 15°C (Fig. [I]), phase sep- 
aration proceeds classically via spinodal decomposition. 
After some time, we obtain two homogeneous phases sep- 
arated by a sharp interface, where 4>i,eq an d 4>2, eq are the 
volume fractions of the upper and lower phase, respec- 
tively, as given by left and right branches of the binodal 
curve in Fig. [T] We proceed with a 'slow' quench in 
this region and observe the spinodal decomposition un- 
der the microscope. At early times, the micrographs in 
Fig. 2a show the typical evolution of a bicontinuous net- 
work with a characteristic size which increases with time, 
indicating a coarsening of the structures formed during 
the spinodal decomposition. After about t ~100 s, the 
dense phase starts to sediment due to the density mis- 
match. Micrographs indeed show an increasingly blurred 
appearance when the focus is moved to the upper part of 
the capillary. From then on the focus is maintained in the 
lower part of the capillary tube. Finally, around £=1000 
s, the dense phase wets the bottom of the capillary and 
spreads to totaly cover the bottom of the capillary tube. 

In parallel USALS experiments were carried out. The 
evolution of the spinodal decomposition process and 
the corresponding characteristic length scale captured in 
Fourier space is shown in Fig. 2b. A peak is observed 
in the scattering pattern at a scattering vector q* that 
corresponds to a characteristic length, £ = 2Tr/q*. With 
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FIG. 2: Microscopy and USALS results obtained during 'clas- 
sical' spinodal decomposition, (a) Micrographs showing the 
time evolution of the spinodal decomposition for a slow shal- 
low quench at 17°C and </>o=0.15. (b) USALS intensity as a 
function of q showing the time evolution of the spinodal de- 
composition for a slow shallow quench at 16°C and </>o=0.11. 
With increasing time the peak moves towards lower (/-values: 
0s (□), lis (A), 17s (o), 34s (V), 101s (o), 118s (<), 151s (>), 
250s (o). (c) Dynamic scaling I{q)q* of the data shown in 
(b). (d) Evolution of the characteristic length obtained dur- 
ing the microscopy and the USALS experiment for the two 
quenches shown in (a) and (c). 



time the peak increases in amplitude and moves towards 
lower scattering vectors, reflecting the coarsening of the 
structures formed during the spinodal decomposition in 
regime (a) already seen in microscopy. At t ^100 s, 
the evolution is perturbed cither by sedimentation of the 
dense phase or by the finite size of the cell. 

Commonly one distinguishes between at least three 
characteristic regimes during the spinodal decomposition 
process [TS]: (i) early, (ii) intermediate (diffusibe) and 
(m) late (flow) stage. The early stage (i) is described by 
linear Cahn-Hillard theory [TT1, HH1 HH] and in essence pre- 
dicts an increase of the amplitude at constant q* until the 
peak spatial concentration fluctuations have reached the 
final coexisting concentrations and domain growth and 
coarsening sets in. This regime is too fast and the struc- 
tural length scale too small to be captured by USALS. 



4 



During the intermediate and late stage, coarsening is ex- 
pected to occur according to a power law £ ~ t a , where 
the scaling exponent is a — 1/3 in the intermediate (ii) 
and a = 1 in the late (Hi) stage of the coarsening regime 
[2TH [2"T] . A detailed calculation for the intermediate stage 
including hydrodynamics reveals that the exponent a is 
between 0.2 and 1.1, depending upon the relative impor- 
tance of hydrodynamic interactions [22 . 

The scattering patterns shown in Fig. 2b are indeed 
typical for spinodal decomposition. This is illustrated 
in Fig. 2c, where the data are shown to follow universal 
scaling by I(q)q* as proposed by Furukawa [33] and de- 
rived by Dhont for the spinodal decomposition of colloids 
in the initial and intermediate stage including hydrody- 
namic interactions [22] . 

A quantitative temporal evolution of £ as obtained 
from the data shown in Fig. 2b through £ = 2ir/q* is 
presented in Fig. 2d. It reflects the three regimes pre- 
viously described: (a) growth (t < 100 s), (b) sedimen- 
tation (100 s< t <10000 s), and (c) wetting-spreading 
(i >10000 s). A comparison with the data from the mi- 
croscopy experiments obtained by calculating the power 
spectrum of the micrographs is also given in Fig. 2d. 
We see that in the intermediate regime, and before the 
onset of perturbations such as sedimentation and finite 
cell size effects, the scattering follows a power law of the 
form £ ~ t a , with a m 1/3, in agreement with theoretical 
predictions and other experiments with colloidal model 
systems, (see [23] and references therein). 



IV. ARRESTED SPINODAL DECOMPOSITION 

Having verified the that the system exhibits all the 
characteristics of classical bicontinuous spinodal decom- 
position for quenches with Tf > 15°, we now follow the 
structural evolution in a sample that is quenched into the 
spinodal region to a final temperature below the arrest 
tie line at 15° C. Here, the sample not only initially un- 
dergoes spinodal decomposition into a bicontinuous net- 
work with protein-rich and protein-poor domains, it also 
exhibits a non-ergodic liquid-solid transition as shown in 
Fig. lb [5]. As these processes turn out to be too fast to 
be monitored by video microscopy and USALS, we focus 
on the final and time-independent structural properties 
of the arrested samples as a function of the parameters 
initial volume fraction (j)Q, quench depth Tt and quench 
rate. 



A. Influence of cj>Q 

We first look at the evolution of the resulting struc- 
ture in samples undergoing an arrested spinodal decom- 
positions as a function of the initial concentration (f>o at 
constant quench temperature Ty=10°C. Fig. 3a presents 
examples of micrographs obtained during 'fast' quenches 
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FIG. 3: Microscopy and USALS results obtained during 
arrested spinodal decomposition at constant quench depth 
T/=10°C and for different (j>o- (a) Micrographs showing ar- 
rested spinodal decomposition for a fast quench at:0o=O.19 
(□),0o=O.15 (o), 00=0.11 (A), 00=0.075 (V). (b) Evolution 
of the characteristic length as a function of <j>o. (o) Results 
for a slow quench obtained during microscopy experiments. 
(•) Results for a slow quench obtained during USALS ex- 
periments. (□) Results for a fast quench obtained during 
microscopy experiments. 



at four different 4>q. While the initial spinodal decom- 
position and subsequent arrest is too fast to be resolved 
and the data collected is time independent, the result- 
ing structures and scattering patterns are very similar 
to the data obtained during 'classical' phase separation. 
We thus proceed with the same analysis to extract £ from 
microscopy and USALS, and the results from these mea- 
surements are summarized in Fig. 3b. Error bars are 
estimated based on 10 successive identical experiments 
with fresh lysozyme samples, and we find£ to vary at 
most within 10% of its mean value. As shown in Fig. 3d, 
the characteristic length depends on the type of quenches, 
but for similar quench rate we observe no variation of £ 
with 4>q. 

The arrest scenario described in |5j as sketched in 
Fig. lc provides an explanation for these findings. In 
the early or linear stage, spinodal decomposition is trig- 
gered by density fluctuations of wavelength A around 4>q 
PH EH [S3J [SB]. The amplitude of the fluctuations in- 
creases with time until the low and high density values 
have reached their equilibrium values given by the corre- 
sponding points on the coexistence curve, i.e., (j>i=4>\^ e q 
and (j)2=4>2,eq- At this point the early stage of the spin- 
odal decomposition ends and the structure starts coars- 
ening. However, the results obtained by [5j indicate 
that in the arrested spinodal decomposition observed for 
lysozyme <f>\ and <f>2 never reach their equilibrium val- 
ues because the dense phase becomes arrested in a glassy 
state at 4>2=4'2,giass < 4>2,eq- Therefore the characteristic 
length we observe corresponds to the the wavelength of 
the density fluctuations. The density fluctuation modes 
grow as exp(Rt) where R is the amplification factor. R 
is proportional to the sum of two antagonistic terms: the 
second derivative of the Helmhotlz energy density of the 
bulk phase Ikf (which is also proportional to the isother- 
mal compressibility) and the energy cost of a concentra- 
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tion gradient k: R ~ ^§ — kc/ . The first term favors 
short wavelengths and it decreases as one approaches the 
spinodal line where it is equal to zero. The second term 
favors long wavelengths to minimize the number of inter- 
faces in the system. The sharpness of the maximum of R 
and the fact that R enters the problem in an exponential 
function imply that the only density fluctuation that is 
relevant is the one for which A maximizes R. Therefore, 
during the early stage of spinodal decomposition it is the 
fastest growing density fluctuation mode of wavelength 

Am ~ \J k/ — that sets the characteristic length of 

the structure: £=A M (I/). In the </> -series, £ does not 
vary with </>o- The likely explanation is that R and thus 
the derivative of the chemical potential and the energy 
cost of a concentration gradient are constant. For those 
values to be constant the local volume fractions 4>i an d 
4>2 must also be constant. This is the case in our exper- 
iments as it is only the total volume of each phase that 
varies with <f>Q\ the total volume of the dense phase in- 
creases and the total volume of the dilute phase decreases 
with 0o- 



B. Influence of T, 



S 



A completely different behavior is observed when we 
look at the evolution of £ for arrested spinodal decompo- 
sition at a constant volume fraction, </>o=0.15 as a func- 
tion of Tf. In Fig. 4a and b we show the micrographs and 
the corresponding £ obtained from a 'fast' quench at four 
different Tf. In Fig. 4b, we also show the evolution of the 
£ as measured by video microscopy and USALS at four 
different Tf for slow quenches. £ obviously depends on 
the quench rate, but for similar quench rates we observe a 
monotonic and pronounced increase of £ with decreasing 
quench depth, i.e., increasing Tf. 

While £ is independent of <j> for fixed quench depth, 
Fig. 4b now clearly demonstrates that it strongly de- 
pends on Tf. We find that £ follows a critical scaling 
of the form £ ~ (T s - Tf) - 55 (Fig.4c) where T s is the 
temperature of the spinodal line. This finding further 
supports our assumption that arrest occurs during the 
early stage of spinodal decomposition. Indeed during the 
early stage of spinodal decomposition the wavelength A M 
of a density fluctuation is proportional to the square root 
of the energy cost of a concentration gradient divided by 
the second derivative of the Hclmhotlz energy density of 
the bulk phase. Am diverges at the spinodal temperature 
because ^| is equal to zero at T s . The increase of £ with 
Tf is qualitatively consistent with the fact that Am di- 
verges at the spinodal temperature T s ((f> = 0.15)=18°C. 
Around the critical volume fraction <j) c , Am follows a crit- 
ical scaling £ ~ (T s - Tff . Although </> =0.15 is slightly 
below the critical volume fraction <f) c ~ 0.17, the expo- 
nent found in our experiments v — 0.55±0.1 is consistent 
with values of the critical exponent v derived from the 
mean field approximation [y = 0.5) and from the renor- 
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FIG. 4: Microscopy and USALS results obtained during ar- 
rested spinodal decomposition at constant 0o=O.15 and for 
different quench depths Tf. (a) Micrographs showing arrested 
spinodal decomposition for a fast quench at: T/ = 15°C (□), 
T/ = 10°C (o), T/=5°C (A), 2>=1°C (V). (b) Evolution of £ 
as a function of Tf. (o) Results for a slow quench obtained 
during microscopy experiments. (•) Results for a slow quench 
obtained during USALS experiments. (□) Results for a fast 
quench obtained during microscopy experiments. The inset 
shows the evolution of the scattered intensity I(q) with tem- 
perature for the four arrested samples measured with USALS, 
where the temperature increases from bottom to top. (c) 
Evolution of the characteristic length as a function of T s — Tf 
where T S =18°C is the temperature of the spinodal at </>o=0.15. 
Lines show linear fit of the data with a slope of 0.6±0.1 for 
(□), 0.55±0.1 for (•) and 0.5±0.1 for (o). The dashed vertical 
line shows the position of the arrest tie line at T/=15°C. 



malization theory (y = 0.63) [151 1251 127] . 

The arrest scenario where the dense phase that forms 
during the early linear stage in the spinodal decomposi- 
tion arrest as its concentration crosses the arrest line is 
further supported by the temperature dependence of I(q) 
(Fig. 4b inset). We observe that the amplitude of the 
peak decreases with decreasing temperature. This im- 
plies that the scattering contrast between the two phases 
is lower. This suggests that the difference between the 
two volume fractions decreases with temperature, which 
is again in agreement with the state diagram in Fig. la 
where the arrest line extends into the unstable region and 
moves to lower volume fractions with decreasing temper- 
atures. 



C. Influence of the quench rate 

The quench rate influences the value of the charac- 
teristic length, with a fast quench leading to smaller £ 
as shown in Fig. 3b and 4b. When quenches are slower, 
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the sample spends more time at temperatures above Tf. 
Therefore density fluctuation modes that correspond to 
these temperatures are activated and contribute to longer 
wavelengths. The combined effects of those different fluc- 
tuation modes results in a larger value of the resulting £. 



V. CONCLUSION 

We have previously shown in ref. [5] that lysozyme 
solutions can undergo an arrested spinodal decomposi- 
tion that results in the formation of solid-like samples. 
The rheological properties of these gels were found to 
exhibit an unusual frequency dependence with two well 
defined plateaus and associated yield stress values. This 
allowed us then to use centrifugation experiments to de- 
termine the local densities of both phases and to precisely 
locate the arrest line close to and within the unstable 
region of the phase diagram. Based on the extension 
of the arrest line in the coexistence region, we assumed 
that arrested spinodal decomposition occurs during the 
early stage when local variations of the density proceed. 
Our measurements now confirm this hypothesis. Using 
a combination of video microscopy and USALS we were 
able to demonstrate the formation of an arrested bicon- 
tinuous network for quenches below the arrest tie line 
at T/=15°C. Moreover, we have been able to show that 
the correlation length in these samples exhibits a tem- 
perature dependence that closely follows critical scaling 
expected for density fluctuations during the early stages 
of spinodal decomposition. 

These findings of course raise new questions. There 
is the apparent inconsistency with the study of Lu ct 
al. [7J |5H] where the authors had proposed that gelation 
for particles with short ranged attraction would happen 
through a spinodal decomposition that proceeds until the 
dense phase undergoes an attractive glass transition at 
(j) g ~ 0.55. One explanation could of course come from 
the fact that while lysozyme exhibits many of the features 
characteristic for colloids with short range attractions, it 
is also known that due to a slightly non-spherical shape 
and patchyness there are some non-centrosymmetric con- 
tributions to the interaction potential that one needs to 
account for in order to quantitatively reproduce the en- 
tire phase diagram |29j . However, an alternative expla- 
nation could come from the way the concentration in the 
arrested phase has been determined by Lu et al. The au- 
thors used a confocal microscope and determined the ac- 
ceptance probability of a test particle randomly inserted 
into their 3D images. This procedure is not unambigu- 
ous for highly ramified systems such as an arrested phase 
formed by attractive particles within a bi-continuous tex- 
ture as seen for lysozyme. Under these conditions the 
method used by Lu et al. would not distinguish between 
space inside and outside the arrested phase, but yield the 
local density of individual gel strands within this phase 
rather then the average concentration in the mesoscopic 
phase. In this context it is also important to note that 



even for the small protein lysozyme with a size of about 
2 nm the extension of the bi-continuous texture as given 
by the correlation length is of the order of a few /mi, i.e. 
about a thousand particles in diameter, and thus much 
larger than the volume in which the (local) concentra- 
tion is determined when using a confocal microscope. It 
is clear that similar measurements as we have done for 
lysozyme would allow to help clarify this issue. 

However, irrespective of the general validity of our find- 
ings with respect to the phase behavior of colloids, our 
results suggest some interesting possibilities. Given the 
strong interest in obtaining gel-like systems with tun- 
able structural and viscoclastic properties for numerous 
applications in areas as diverse as materials and food sci- 
ence, arrested spinodal decomposition where the solidifi- 
cation occurs in the early stage offers means to create gels 
with a full control over the resulting correlation length of 
the bi-continuous network. This work also supports and 
gives new perspectives to the recent work of E. Dufrcsnc 
where they have hypothesize that multiple lineages of 
birds have convergently evolved to exploit phase sepa- 
ration and kinetic arrest to self-assemble spongy color- 
producing nano-structures in feather barbs [30] . More- 
over, as we already know from our previous rheological 
measurements that the arrested networks have a material 
response function that is characterized by two distinct 
elastic moduli that we believe to be coupled to the two 
characteristic length scales £ and protein diameter [5], 
the obtainable variability in £ may result in large vari- 
ations of the viscoelastic properties of these soft solids, 
and thus offer interesting possibilities for creating pro- 
tein gels with tailored texture. We thus currently work 
on a full characterization of the viscoelastic properties 
of systems that result from different quenches and thus 
different structural properties. 
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